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Abstract 



We discuss the properties of an ideal relativistic gas of events possessing 
Bose-Einstein statistics. We find that the mass spectrum of such a system is 
bounded by fi < m < 2M/fj,x, where fi is the usual chemical potential, M is 
an intrinsic dimensional scale parameter for the motion of an event in space- 
time, and hk is an additional mass potential of the ensemble. For the system 
including both particles and antiparticles, with nonzero chemical potential //, 
the mass spectrum is shown to be bounded by < m < 2M/ fix, and a special 
type of high-temperature Bose-Einstein condensation can occur. We study this 
Bose-Einstein condensation, and show that it corresponds to a phase transition 
from the sector of continuous relativistic mass distributions to a sector in which 
the boson mass distribution becomes sharp at a definite mass M/fix- This 
phenomenon provides a mechanism for the mass distribution of the particles to 
be sharp at some definite value. 

Key words: special relativity, relativistic Bose-Einstein condensation, mass distribu- 
tion, mass shell 
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1 Introduction 

There have been a number of papers in the past |], 0, |3|, U > which discuss the properties 
of an ideal relativistic Bose gas with nonzero chemical potential fi. Particular attention 
has been given to the behavior of the Bose-Einstein condensation and the nature of 
the phase transition in d space dimensions @, |^] . The basic work was done many years 
ago by Jiittner ||, Glaser J7|, and more recently by Landsberg and Dunning-Davies 
|| and Nieto 0. These works were all done in the framework of the usual on-shell 
relativistic statistical mechanics. 

To describe an ideal Bose gas in the grand canonical ensemble, the usual expression 
for the number of bosons N in relativistic statistical mechanics is 



where V is the system's three-volume, Ek = Vk 2 + m 2 and T is the absolute tem- 
perature (we use the system of units in which ft = c = /eg = 1; we also use the 
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metric g^ u = (—,+,+,+)), and one must require that // < m in order to ensure a 
positive-definite value for n&, the number of bosons with momentum k. Here N is 
assumed to be a conserved quantity, so that it makes sense to talk of a box of N 
bosons. This can no longer be true once T ~ m [10]; at such temperatures quantum 
field theory requires consideration of particle-antiparticle pair production. If N is the 
number of antiparticles, then N and N by themselves are not conserved but N — N 
is. Therefore, the high-temperature limit of (1.1) is not relevant in realistic physical 
systems. 

The introduction of antiparticles into the theory in a systematic way was made 
by Haber and Weldon [It], [TT]] . They considered an ideal Bose gas with a con- 
served quantum number (referred to as "charge") Q, which corresponds to a quan- 
tum mechanical particle number operator commuting with the Hamiltonian HQ All 
thermodynamic quantities may be then obtained from the grand partition function 
Tr {exp [—(H — fiQ)/T]} considered as a function of T, V, and fi [|12|| . The formula 
for the conserved net charge, which replaces (1.1), readsQ [10] 



q = vj: 



i 



e {E k -y)/T _ I e (E k +n)/T _ I 
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In such a formulation a boson-antiboson system is described by only one chemical 
potential //; the sign of p, indicates whether particles outnumber antiparticles or vice 
versa. The requirement that both and be positive definite leads to the important 
relation 

\n\ < m. (1.3) 

The sum over k in (1.2) can be replaced by an integral, so that the charge density 
p = Q/V becomes 
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e {E k -ti)/T _ I 
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which is an implicit formula for \x as a function of p and T, and in the region T » m 
reduces to 

P*>£. (1-5) 

For T above some critical temperature T c , one can always find a fi (|//| < m) such 
that (1.4) holds. Below T c no such /i can be found, and (1.4) should be interpreted 
as the charge density of the excited states: p — po, where po is the charge density 
of the ground state [10] (with k = 0; clearly, this state is given with zero weight in 

1 In the manifestly covariant theory which we shall use in our study, this charge is naturally 
associated with particles and antiparticles which are distinguished by the off-shell structure, as in 



quantum field theory [10 . 

2 The standard recipe according to which all additive thermodynamic quantities are reversed for 
antiparticles is used. 
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the integral (1.4)). The critical temperature T c at which Bose-Einstein condensation 
occurs corresponds to p = ±m (depending on the sign of p.) Thus, one sets |p| = m 
in (1.4) and obtains, via (1.5) (provided that \p\ » m 3 ), 



I3\p\ 
m 



;i.6) 



Below T c , (1.4) is an equation for p — p , so that the charge density in the ground 
state is 

Po =p[l-(T/T c ) 2 ]. (1.7) 

It follows from Eq. (1.6) that any ideal Bose gas will condense at a relativistic 
temperature (T c >> m), provided that |p| >> m 3 . 

Recently an analogous phenomenon has been studied in relativistic quantum field 
theory (jy], [13|, [14], [15|] . For relativistic fields Bose-Einstein condensation occurs at high 



temperatures and can be interpreted in terms of a spontaneous symmetry breaking 
0. 

In this paper we shall use a manifestly covariant form of statistical mechanics 
which has more general structure than the standard forms of relativistic statistical 
mechanics, but which reduces to those theories in a certain limit, to be described 
precisely below. In fact, it is one of the principle aims of this work to provide a 
mechanism for which this limit can be realized on a statistical level. The results that 
we obtain are different from those of the standard theories at high temperatures. 
These theories, which are characterized classically by mass-shell constraints, and the 
use, in quantum field theory, of fields which are constructed on the basis of on- 
mass-shell free fields, are associated with the statistical treatment of world lines and 
hence, considerable coherence (in terms of the macroscopic structure of whole world 
lines as the elementary objects of the theory) is implied. In nonrelativistic statistical 
mechanics, the elementary objects of the theory are points. The relativistic analog of 
this essentially structureless foundation for a statistical theory is the set of points in 
spacetime, i.e., the so-called events, not the world lines (Currie, Jordan and Sudarshan 
|i~6| have discussed the difficulty of constructing a relativistic mechanics on the basis 
of world lines). 

The mass of particles in a mechanical theory of events is necessarily a dynamical 
variable, since the classical phase space of the relativistic set of events consists of 
the spacetime and energy-momentum coordinates {q«,ti; Pi,Ei}, with no a priori 
constraint on the relation between the p« and the E i: and hence such theories are 
"off-shell". It is well known from the work of Newton and Wigner jTT]] that on- 
shell relativistic quantum theories such as those governed by Klein-Gordon or Dirac 
type equations do not provide local descriptions (the wave functions corresponding 
to localized particles are spread out); for such theories the notion of ensembles over 
local initial conditions is difficult to formulate. The off-shell theory that we shall use 



here is, however, precisely local in both its first and second quantized forms [18, 19 
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The phenomenologial predictions of on-shell theories, furthermore, provide equa- 
tions of state which appear to be too rigid. Shuryak has obtained equations 
of state which are more realistic by taking into account the spectrum of mass as 
seen in the resonance spectrum of strongly interacting matter. We have shown [[H|] 
that Shuryak's "realistic" equation of state follows in a natural way from the mass 
distribution functions of the off-shell theory. 

We finally remark that the standard formulations of quantum relativistic sta- 
tistical mechanics, and quantum field theory at finite temperature, lack manifest 
covariance on a fundamental level. As for nonrelativistic statistical mechanics, the 
partition function is described by the Hamiltonian, which is not an invariant ob- 
ject, and hence thermodynamic mean values do not have tensor properties. [One 
could consider the invariant p^n^ in place of the Hamiltonian [22], where is a unit 



four-vector; this construction (supplemented by a spacelike vector othogonal to n^) 
implies an induced representation for spacetime. The quantity that takes the place 
of the parameter t is then x^n^. This construction is closely related to the problem 



pointed out by Currie, Jordan and Sudarshan ||16fl , for which different world lines 
are predicted dynamically by the change in the form of the effective Hamiltonian 
in different frames.] Since the form of such a theory is not constrained by covari- 
ance requirments, its dynamical structure and predictions may be different than for 
a theory which satisifies these requirements. For example, the canonical distribution 



of Pauli [E3| for the free Boltzmann gas has a high temperature limit in which the 
energy is given by 3ksT, which does not correspond to any known equipartition rule, 
but for the corresponding distribution for the manifestly covariant theory, the limit 
is 2k bT, corresponding to \k B T for each of the four relativistic degrees of freedom. 
For the quantum field theories at finite temperature, the path integral formulation 
24j replaces the Hamiltonian in the canonical exponent by the Lagrangian due to 



the infinite product of factors (4>\tt) (transition matrix element of the canonical field 
and its conjugate required to give a Weyl ordered Hamiltonian its numerical value). 
However, it is the t variable which is analytically continued to construct the finite 
temperature canonical ensemble, completely removing the covariance of the theoret- 
ical framework. One may argue that some frame has to be chosen for the statistical 
theory to be developed, and perhaps even for temperature to have a meaning, but 
as we have remarked above, the requirement of relativistic covariance has dynamical 
consequences (note that the model Lagrangians used in the non-covariant formula- 
tions are established with the criterion of relativistic covariance in mind), and we 
argue that the choice of a frame, if necessary for some physical reason, such as the 
definition and measurement of temperature, should be made in the framework of a 
manifestly covariant structure. 

We consider, in this paper, a relativistic Bose gas within the framework of a 



manifestly covariant relativistic statistical mechanics ER 26, E7]]. We obtain the 



expressions for characteristic thermodynamic quantities and show that they coincide 
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quantitatively, in the narrow mass-width approximation, with those of the relativistic 
on-shell theory, except for the value of the average energy (which differs by a factor 
2/3, as remarked above). We introduce antiparticles and discuss the high temperature 
Bose-Einstein condensation in such a particle-antiparticle system. We show that it 
corresponds to a phase transition to a high-temperature form of the usual on-shell 
relativistic kinetic theory. In the following, we briefly review the manifestly covariant 
mechanics and quantum mechanics which forms the basis of our study of relativistic 
statistical mechanics. 

In the framework of a manifestly covariant relativistic statistical mechanics, the 
dynamical evolution of a system of N particles, for the classical case, is governed by 
equations of motion that are of the form of Hamilton equations for the motion of N 
events which generate the space-time trajectories (particle world lines) as functions 



of a continuous Poincare-invariant parameter r, called the "historical time" [p8|, |29 
These events are characterized by their positions g M = (t, q) and energy-momenta 
p^ = (E, p) in an 8iV-dimensional phase-space. For the quantum case, the system 
is characterized by the wave function ip T (qi, <?2> • • • , Qn) £ L 2 (R m ), with the measure 
d A qid i q2 • ■ ■ d A q^ = d m q, (g, = gf ; /x = 0, 1, 2, 3; i = 1,2,..., AT), describing the 
distribution of events, which evolves with a generalized Schrodinger equation ||29|| . The 



collection of events (called "concatenation" j3(J) along each world line corresponds 
to a particle, and hence, the evolution of the state of the iV-event system describes, 
a posteriori, the history in space and time of an iV-particle system. 



For a system of N interacting events (and hence, particles) one takes |29] 



K = Y,% L + V(quq2,...,q N ), (1.8) 

i 

where M is a given fixed parameter (an intrinsic property of the particles), with the 
dimension of mass, taken to be the same for all the particles of the system. The 
Hamilton equations are 

dqj = dK_ = pf_ 
dr dp ifi M' 

drf OK dV 
dr dq ilM dq ifl ' 
In the quantum theory, the generalized Schrodinger equation 



;i.9) 



d 

i-Q-^AQu 92, • • • , ?tv) = K^ T (q u q 2 ,..., q N ) (1.10) 

describes the evolution of the iV-body wave function ip T (qi, q 2 , ■ ■ ■ , gjv)- To illustrate 
the meaning of this wave function, consider the case of a single free event. In this 
case (1.10) has the formal solution 

Mq) = (e-^oXg) (i- 11 ) 
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for the evolution of the free wave packet. Let us represent ip T {q) by its Fourier 
transform, in the energy-momentum space: 

Mo) = 777^? / d 4 pe-^ T e^Mp), (1.12) 



(2tt)' 

where p 2 = p^p^, p-q = p^q^, and ipo(p) corresponds to the initial state. Applying the 
Ehrenfest arguments of stationary phase to obtain the principal contribution to ip T (q) 
for a wave packet at p%, one finds (p£ is the peak value in the distribution ipo(p)) 

€ * §r, (1.13) 

consistent with the classical equations (1.9). Therefore, the central peak of the wave 
packet moves along the classical trajectory of an event, i.e., the classical world line. 

In the case that p° c = E c < 0, we see, as in Stueckelberg's classical example ||28|| , 
that 

dt c E c 

— — — -r-z < 0. 

dr M 

It has been shown f3(| in the analysis of an evolution operator with minimal electro- 
magnetic interaction, of the form 

(p - eA{q)f 
2M 

that the CPT-conjugate wave function is given by 

Vv CPT (t,q)=VvH,-q), (1.14) 
with e — > —e. For the free wave packet, one has 

^ PT (q) = T^p/ d'pe-^e-^Mp)- (1-15) 
The Ehrenfest motion in this case is 

if E c < 0, we see that the motion of the event in the CPT -COHJUgclt6 stclt6 is in the 
positive direction of time, i.e., 



dt c E c \E C \ 



dr M M 



;i.i6) 



and one obtains the representation of a positive energy generic event with the opposite 
sign of charge, i.e., the antiparticle. 
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It is clear from the form of (1.10) that one can construct relativistic transport the- 
ory in a form analogous to that of the nonrelativistic theory; a relativistic Boltzmann 
equation and its consequences, for example, was studied in ref. [^Bj . 



As a simple example of the implications of the classical dynamical equations (1.9), 
consider the problem of a relativistic particle in a uniform external "gravitational" 
field, with evolution function 

K = p I y_ + M / L17 ) 

(the external potential breaks the invariance of the evolution function, but that will 
not affect the illustrative value of the example) with initial conditions t(0) = 0, t(0) = 
a, z(0) — h, i(0) = 0, resulting in the solution 

1 2 , 

z = ~-j9 T t = ar + t , 

E = Mc 2 a, p z = -Mgr. (1.18) 

The invariant variable r replaces t in describing the dynamical evolution of the system. 
The generator of the motion 

v 2 - E 2 lc 2 1 
K = — — — + mgz = -Mc 2 a 2 = const, (1.19) 

as required. The total energy of the particle in this case, including both increase of 
momentum and decrease of dynamical mass, is constant also. The effective particle 
mass m is given by 



m 



y { E/cy- P l = MaJl- 9 ^. (1.20) 



Expanding this out in the nonrelativistic limit c —>■ oo, one obtains (with r 2 = 
2(h - z)/g) 

m = Ma- — z), (1.21) 

ac 1 

and we recognize Mgih — z)/c 2 as the mass shift induced by the potential term. 
The factor a arises due to the choice of initial conditions, i.e., for r = 0, m = Ma, 
and not M (for r sufficiently large, under this unbounded potential, the quantity in 
the square root could become negative, and the particle could become tachyonic). 
Note that it is the mass of the particle which carries dynamical information (the 
total energy is constant, but the mass is "redshifted" by the potential) and that has 
the correspondence with nonrelativistic energy, through the mass-energy equivalence, 
that we observe in the laboratory. This point is discussed in more detail in, for 
example, refs. [31] and 
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2 Ideal relativistic Bose gas without antiparticles 



To describe an ideal gas of events obeying Bose-Einstein statistics in the grand canon- 
ical ensemble, we use the expression for the number of events found in |25 



iv = ^E^ = ^ (4) E' 



i 



k>* e 



(2.1) 



where is the system's four-volume and m 2 = —A; 2 = —k^k^; jix is an additional 



mass potential [p5 |, which arises in the grand canonical ensemble as the derivative 
of the free energy with respect to the value of the dynamical evolution function K, 
interpreted as the invariant mass of the system. In the kinetic theory |25], enters 
as a Lagrange multiplier for the equilibrium distribution for K, as \i is for N, and 1/T 
for E. We shall see, in the following, how fix plays a fundamental role in determining 
the structure of the mass distribution. In order to simplify subsequent considerations, 
we shall take it to be a fixed parameter. 

To ensure a positive-definite value for n^n, the number density of bosons with 
four-momentum k^, we require that 



m 

m — a — a K > 0. 

The discriminant for the l.h.s. of the inequality must be nonnegative, i.e., 

M 



< 



2fi K 



For such fi, (2.2) has the solution 



mi 




m 2 . 



For small fifix/M, the region (2.4) may be approximated by 



jit < m < 



2M 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



One sees that [lk determines an upper bound of the mass spectrum, in addition to 
the usual lower bound m > /i. In fact, small fix admits a very large range of off-shell 
mass, and hence can be associated with the presence of strong interactions |33] . 



Replacing the sum over k^ (2.1) by an integral, one obtains for the density of 
events per unit space-time volume n = N/V {i) §3, 



n 



An 3 



1112 



dm 



dp 



m 3 sinh j3 



(2.6) 
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where mi and iri2 are defined in Eq. (2.4), and we have used the parametrization [f2£j 

p° = m cosh j3, 

p 1 = m sinh (3 sin # cos </>, 

p 2 = m sinh /3 sin # sin 0, 

p 3 = m sinh /3 cos 0, 

< < 7T, < <p < 2tt, -oo < {3 < oo. 
In this paper we shall restrict ourselves to the case of high temperature alone: 

M , . 

T » — . (2.7) 

It is then possible to use, for simplicity, the Maxwell-Boltzmann form for the inte- 
grand, and to rewrite (2.6) in the form 

n = 



Arc 3 



/ m 3 dm sinh 2 /?rf/?e- mcosh/3/T e^ m /2MT , (2.8) 

Jm\ J —oo 



which reduces, upon integrating out /?, to P7 



471-3 



r dm m 2 K x ( e M— 2 /2MT ) (2 9) 



where K v {z) is the Bessel function of the third kind (imaginary argument). Since 
/i < m < m 2 < 2M//ik, 

/W < fi K (2M/fi K ) 2 = _2M j f2 im 

2MT ~ 2MT T/i X ' 1 ' J 

in view of (2.7), and also 

a m 2M 

Therefore, one can neglect the exponentials in Eq. (2.9), and for Ki(m/T) use the 
asymptotic formula [35 

Kv{z)~lr(v)(Zy V , z«l. (2.12) 

Then, we obtain 



2 



n ~£r^ m= *L "YM-lgL. (2.13) 



4tt 3 Jmi 2vr 3 \(i K J V M 
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^From this equation, one can identify the high-temperature mass distribution for the 
system we are studying, so that now 



, (1 C>dmm e+1 _ 2 m e 2 +2 -m[ +2 

" m) ~ Qdmm ~ i + 2 m 2 2 -ml ' (2-J1) 



W = 77 1 - 7T7T > ( 2 - 15 ) 



In particular, 

(771.) = 

V 2M J 

Extracting the joint distribution for /? and m from (2.8) in the same way, we also 
obtain the average values of the energy and the square of the energy for high T. The 
average energy is given by 

, , . . J™ 2 m 4 dmsmh. 2 pcoshf3d(3e- mcoshl3 / T 

(E) = (mcosh/3) - Cm3dmsinh 2 We _ mMBh , /T • (2-17) 

Integrating out /3, one finds 

(F) ^ I C dm m^[K 3 (m/T) - K^m/T)] 

[ ' AT SZZdmmtK^m/T) ' 1 ' 

It is seen, with the help of (2.12), that it is possible to neglect K\ in comparison with 
K 3 in the numerator of (2.18) and obtain, via (2.12), 

, E , 1 f%dmm*K 3 (m/T) 

^ ; ~ ATJ^dmm^K^m/T) ' 1 J 

in agreement with refs. |25, ^6], pTfl . Similarly, one obtains 

J™ 2 m 5 dm sinh 2 (3 cosh 2 f3df3e- mcosh/3 / T 



(E ) = (m 2 cosh 2 /3) 



/™ 2 m 3 dm sinh 2 f3d(3e~ 



m cosh /3/T 



= C dm[m%(m/T) + 3Tm 3 iT 2 (m/T)] ^ ^ ffi dmm*K 2 (m/T) ^ ^ 
f™ 2 dm m 2 Ki(m/T) J™ 2 dm m 2 Ki(m/T) 

(2.20) 

Let us assume, as is generally done, that the average (p^p") has the form 

(p/y) = aw"/ + bg^, (2.21) 

where m m = (1, 0) in the local rest frame. The values of a and b can then be calculated 
as follows: for \x = v = one has ((p ) 2 ) = a — b, while contraction of (2.21) with g^ u 
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gives —g^ u (p^Pu) = a — 4b. The use of the expressions (2.20) for ((p ) 2 ) = (E 2 ), and 
(2.16) for —g^ipfiPv) = {m 2 ) yields 

J a -b = 6T 2 , 
i o-46 = 2(^) 2 (1-/WM), 



so that 



2 ( M\ 2 f fifi K \ 



6 - 2r -aUJ I 1 " m J- < 2 ' 23 > 

For T >> M/fix, it is possible to take a = 8T 2 , 6 = 2T 2 , and obtain, therefore, 

(pY) = 8T 2 uV + 2T 2 g flu . (2.24) 

To find the expressions for the pressure and energy density in our ensemble, we 
study the particle energy- momentum tensor defined by the relation 

T*"(q) = £ / drj^-5\q - ft (r)), (2.25) 



in which M/hk is the value around which the mass of the bosons making up the 
ensemble is distributed, i.e., it corresponds to the limiting mass-shell value when the 
inequality (2.3) becomes equality. Upon integrating over a small space-time volume 
AV and taking the ensemble average, (2.25) reduces to 



(Tn = T^ n ^P v ). (2.26) 

In this formula T&y is the average passage interval in r for the events which pass 
through the small (typical) four-volume AV^ in the neighborhood of the _R 4 -point. The 
four- volume AV is the smallest that can be considered a macrovolume in representing 
the ensemble. Using the standard expression 

(7W) = (p + P)u»u v + pg^, (2.27) 

where p and p are the particle pressure and energy density, respectively, we obtain 



p ^ pM = ^f^Jl- 2 -tBL T ', p = 3p. (2.28) 

To interpret these results we calculate the particle number density per unit three- 
volume. The particle four-current is given by the formula p6l 



■ n "^^S ,h i£-/ (q - q - (T)l (2 - 29) 
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which upon integrating over a small space-time volume and taking the average reduces 
to 

(J") = t^-^>; (2-30) 



then 



so that 



No = (J°) = i^n(E), (2.31) 

M/fJL K 



and we recover the ideal gas law 

p = NqT. (2.33) 

Since, in view of (2.4), 

2M [ 2^ 

1/1 n— = ^ m 

is a width of the mass distribution around the value M/fi K , Eqs. (2. 28), (2. 32) can 
be rewritten as 

T AV Am 4 

N ° = o 3 T - 2 - 34 



In ref. [36] we obtained the formulas for thermodynamic variables, under the assump- 
tion of narrow mass width, which depend on T A yAm as well; the requirement that 
these results coincide with those of the usual on-shell theories implies the relation^ 

T AV Am = 2tt. (2.35) 

One can understand this relation, up to a numerical factor, in terms of the uncertainty 
principle (rigorous in the L 2 (R 4 ) quantum theory) AE ■ At ~ 1/2. Since the time 
interval for the particle to pass the volume AV (this smallest macroscopic volume 
is bounded from below by the size of the wave packets) At = E/M At, and the 
dispersion of E due to the mass distribution is AE ~ mAm/E, one obtains a lower 
bound for T AV Am of order unity. 

Thus, with (2.35) holding, the formulas (2.34) reduce to 

V = ~2, P = 3p, (2.36) 
N = ^, (2.37) 



3 In c.g.s. units, this relation has a factor h/c 2 on the right hand side. 
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which are the standard expressions for high temperature |J7| . The formulas for char- 
acteristic thermodynamic quantities and the equation of state for a relativistic gas of 
off-shell events have the same form as those of the relativistic gas of on-shell particles. 
They coincide with them (under the condition (2.35)) in the narrow mass shell limit, 
except for the expression for the average energy which takes the value 2T in the rela- 
tivistic gas of events, in contrast to 3T, as for the high-temperature limit of the usual 



theory ||23|| . Experimental measurement of average energy at high temperature can, 
therefore, affirm (or negate) the validity of the off-shell theory. There seems to be 
no empirical evidence which distinguishes between these results at the present time. 
The quantity o = Mqc 2 /ksT, a parameter which distinguishes the relativistic from 
the nonrelativistic regime (see, e.g., [p5fl ) is very large for M of the order of the pion 
mass, at ordinary temperatures; the ultrarelativistic limit corresponding to a small 
becomes a reasonable approximation for T ~ 10 12 K. 



3 Antiparticles and condensation 



The introduction of antiparticles into the theory as the CPT conjugate of negative 
energy events leads, by application of the arguments of Haber and Weldon 
Actor 
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38 



, or 
We 



to a change in sign of /i in the distribution function for antiparticles 
therefore write down the following relation which represents the analog of the formula 
(1.2): Q 

1 1 



N = v {4) 



(3.1) 



With respect to the determination of the sign of the second term, let us consider a 
space-time picture in which we have many world lines, generated by events moving 
monotonically in the positive t direction. The addition of a particle-antiparticle pair 
which annihilates corresponds to the addition of a world line which is generated by 
an event initially moving in the positive direction of time to some upper bound, to, 
where annihilation takes place, and returning in the negative direction of time. At 
times later than to the total particle number is unaffected. At times earlier than 
to, a particle and antiparticle are added to the total paticle number. Since, as also 
assumed by Haber and Weldon [TO] , the total particle number is a conserved quantity, 
the antiparticle trajectory must be counted with a sign oppostie to that of the particle 
trajectory. The second term in (3.1), counting antiparticles, must therefore carry a 

4 As for the nonrelativistic theory, the "free" distribution functions describe quasiparticles in a 
form which takes interactions into account entering through the chemical potential. By definition, 
good quasiparticles are not frequently emitted or absorbed; we therefore consider the (quasi-) par- 
ticles and antiparticles as two species. Since the particle number is determined by the derivative of 
the free energy with respect to the chemical potential, \i must change sign for the antiparticles flo|| . 
Similarly, the average mass (squared) is obtained by the derivative with respect to [Ik |p5[ ; since 
the mass (squared) of the antiparticle is positive, [Lk does not change sign. 
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negative sign. We require that both terms in Eq. (3.1) be positive definite. In 
this way we obtain the two quadratic inequalities, 



2 

m 

m ~ fI ~ fiK 2M ~ °' 

o 

'ffl 

m + n- Hk— > 0, (3.2) 

which give the following relation representing the nonnegativeness of the correspond- 
ing discriminants: 

MM , x 

<u< . 3.3 

2fi K ~ * ~ 2^ K K J 

It then follows that we must consider the intersection of the ranges of validity of the 
two inequalities (3.2). Indeed, if each inequality is treated separately, there would be 
some values of m for which one and not another would be physically acceptable. One 
finds the bounds of this intersection region by solving these inequalities, and obtains^ 



M 2\ii\ii k \ M I 2\h\hk \ 

which for small |/i|/i^/M reduces, as in the no-antiparticle case (2.5), to 

< m < . (3.5) 

Replacing the summation in (3.1) by integration, we obtain a formula for the 
number density: 



47T 3 



— / m 3 dm / sinh 2 (3d(3 

7T Jmi J — oo 



1 1 



e (mcosh P-fi-fi K ^)/T _ j e (moaab.p+ft-n K ^)/T _ ]_ 

(3-6) 

where mi and m 2 are defined in Eq. (3.4), which for large T reduces, as above, to 



e*i — e 

n = 



471-3 



r nio / 7Tl\ 

-T / dm m 2 K 1 I — ) 

J mi V T ) 



Now, using the estimates (2. 10), (2. 11), and sinh(/i/T) = ji/T for /i/T << 1, we 
obtain (in place of (2.13)) the net event charge 



5 This is actually the solution of one of the inequalities (3.2) (the most restrictive), depending on 
the sign of fj,. 
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The pressure and energy density are obtained by the sum particle and antiparticle 
contributions (proportional to exp(±///T), with the number density (3.7). To second 
order in (/i/T) 2 , one finds 

P = 2 P(H); 
P = 2p(|//|), 

where p(fi) and are given by (2.28) with \i replaced by On the other hand, 
from (2.31) and (3.7), one finds 



Tav M / 2\u\uk m 9 
No = 2^—1/1 " Jl jr L pT 2 , 3.8 

where the factor of 2/i/T, as compared to (2.32), arises from the difference between 
the factors exp(±/x/T) (the sign of \x indicates whether particles or antiparticles 
predominate). One then obtains the following expressions for the Bose gas including 
both particles and antiparticlesf] (here, Am is not necessarily small): 



T AV Am 2T 4 

P = o 5"' ^ = 3 - 9 

27T 7T 

Ta\/ Am 2T 2 

iV = — (3.10) 

Z7T 7T Z 



We now wish to show that the dynamical properties of the current, which follow 
from the relativistic canonical equations of motion, are consistent with the thermo- 
dynamic relation 

N 

where N is the number of bosons in a three-dimensional box of volume V. Since the 
event number density n is, by definition, 

N N 

n 



VAt' 

where At is the (average) extent of the ensemble along the g°-axis (as in our discussion 
after (2.35)), one has 

N = nAt. (3.12) 

6 If we did not neglect indistinguishability of bosons at high temperature, we would obtain, instead 
of (2.37) §fj, N = 2^Zi 3 ( e ' i / T ), where Li v (z) = J2T=i z * I sV is the polylogarithm |9), so that, 
for the system including both particles and antiparticles, No = ^■[Lis^e^ 7 ') — Li^{e~ ll / T )\. It then 
follows from the properties of the polylogarithms J|§] that, for x = \fi\/T « 1, Li s (e x )~ Li 3 (e~ x ) = 
\x, so that, we would obtain, instead of (3.10), N = [iT 2 /3, which coincides with Haber and 
Weldon's Eq. (1.5). 
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The equation of motion (1.9) for q° (with M/fiK, the central value of the mass dis- 
tribution, instead of M, which corresponds to a change of scale parameter in the 
expression (1.8) for the generalized Hamiltonian K), 

dq° _ P° 
dr Mj hk ' 

upon averaging over the whole ensemble, reduces to 

At (E) 

T~ = xb^> 3 ' 13 
T A v M/hk 

where T^v is the average passage interval in r used in previous consideration. Then, 
in view of (3. 12), (3. 13), one obtains the equation (2.31). 



3.1 Relativistic Bose-Einstein condensation 

Since in the particle-antiparticle case, N re \ = N — N, where N and N are the numbers 
of particles and antiparticles, respectively, is a conserved quantity, according to the 
arguments of Haber and Weldon [HJ pointed out in Section 1, and our discussion 
above, N = N ie \/V is also a conserved quantity, so that it makes sense to talk of 
|iV rel | bosons in a spatial box of the volume V. Therefore, in Eq. (3.10) N is a 
conserved quantity, so that, the dependence of /i on temperature is defined by (we 
assume that N is continuous at the phase transition) 

27T 7T 2 iVo 

V=r — X T^pr- 3 - 14 

1 Ay Am 21 1 

For T above some critical temperature, one can always find \i satisfying (3.3) 
such that the relation (3.14) holds; no such \l can be found for T below the critical 
temperature. The value of the critical temperature is defined by putting |/i| = M/2/i^ 
in (3.14). In the narrow mass-shell limit, inserting (2.35), one obtains 



|JVb| 



\ Mj ILK 



(3.15) 



For = M/2fiK, the width of the mass distribution is zero, in view of (3.4), and 
hence the ensemble approaches a distribution sharply peaked at the mass-shell value 
M/fLK- The fluctuations 5m = y (m 2 ) — (m) 2 also vanish. Indeed, as follows from 
(2.15), (2. 16) with /i replaced by \/l\, and (3. 14), (3.15), 
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so that, at T = T c , 5m = 0. It follows from (3.16) that for T in the vicinity of T c 

(T > T c ), 




as for a second order phase transition, for which fluctuations go to zero smoothly. 

We note that Eqs. (2. 36), (2. 37) do not contain explicit dependence on the chemical 
potential, and hence no phase transition is induced. In fact, at lower temperature (or 
small fix) one or the other of the particle or antiparticle distribution dominates, and 



one returns to the case of the high-temperature strongly interacting gas |40| . The 
remaining phase transition is the usual low-temperature Bose-Einstein condensation 
discussed in the textbooks. 

One sees, with the help of (3.4), that the expression for n (3.7) can be rewritten 

as 

1 M A t \ 

n = — - — AmuT; 3.18 

27T d H K 

since at T = T c , Am = 0, it follows that n = at all temperatures below T c . 
Therefore, the behaviour of an ultrarelativistic Bose gas including both particles and 
antiparticles, which is governed by the relation (3.14), can be thought of as a special 
type of Bose-Einstein condensation to a ground state with p^p^ = —(M/hk) 2 (this 
ground state occurs with zero weight in the integral (3.6)). In such a formulation, 
every state with temperature T > T c , given by Eq. (3.6), should be considered as 
an off-shell excitation of the on-shell ground state. At T = T c , all such excitations 
freeze out and the distribution becomes strongly peaked at a definite mass, i.e., the 
system undergoes a phase transition to the on-shell sector. Note that, for n = 0, Eq. 
(3.12) gives At = oo. Then, since (E) ~ T, one obtains from (3.13) that T&y = oo 
(this relation can be also obtained from (2.35) for Am = 0), which means that in the 
mean, all the events become particles. 

As the distribution function enters the on-shell phase at T = T c , the underlying 
off-shell theory describes fluctuations around the sharp mean mass. This phenomenon 
provides a mechanism, based on equilibrium statistical mechanics, for understanding 
how the general off-shell theory is constrained to the neighborhood of a sharp universal 
mass shell for each particle type. At temperatures below T c , the results of the theory 
for the main thermodynamic quantities coincide with those of the usual on-shell 
theories. 

In order that our considerations be valid, the relation T c » M/hk must hold; 
this relation reduces, with (3.15), to 



1 / M N 



»- — . (3.19) 



7T 2 \HK 



For M/fi K ~ ~ 140 MeV, this inequality yields iVo >> 3 • 10 5 MeV 3 . Taking 
iVo ~ 5 • 10 6 MeV 3 , which corresponds to temperature ~ 350 MeV, in view of (2.37), 
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one gets T c ~ 550 MeV ~ Am^. 

If fix is very small, it is difficult to satisfy (3.19) and the possibility of such a 
phase transition may disappear. This case corresponds, as noted above, to that of 
strong interactions and is discussed in succeeding paper p0 . 



4 Concluding remarks 



We have considered the ideal relativistic Bose gas within the framework of a manifestly 
covariant relativistic statistical mechanics, taking account of antiparticles. We have 
shown that in such a particle-antiparticle system at some critical temperature T c a 
special type of relativistic Bose-Einstein condensation sets in, which corresponds to 
a phase transition from the sector of relativistic mass distributions to a sector in 
which the boson mass distribution peaks at a definite mass. The results which can be 
computed from the latter coincide with those obtained in a high-temperature limit of 
the usual on-shell relativistic theory. 

The relativistic Bose-Einstein condensation in particle-antiparticle system consid- 
ered in the present paper can represent (as for the Galilean limit c — > oo ||36|| ) a 
possible mechanism of acquiring a given sharp mass distribution by the particles of 
the system, as a phase transition between the corresponding sectors of the theory. 
Since this phase transition can occur at an ultrarelativistic temperature, it might be 
relevant to cosmological models. The relativistic Bose-Einstein condensation consid- 
ered in the present paper may also have properties which could be useful in the study 
of relativistic boson stars [41]. These and the other aspects of the theory are now 
under further investigation. 

The extension and generalization of Bose-Einstein condensation to curved space- 
times and space-times with boundaries, for which the work reported here may have 
constructive application, has also been the subject of much study. The non-relativistic 
Bose gas in the Einstein static universe was treated in ref. 0. The generalization to 
relativistic scalar fields was given in refs. E2L E3J. The extension to higher dimensional 
spheres was given in ref. [44|. Bose-Einstein condensation on hyperbolic manifolds 

m 



p[5| , and in the Taub universe 



has also been considered. More recently, by 
calculating the high-temperature expansion of the thermodynamic potential when the 
boundaries are present, Kirsten ||7[] examined Bose-Einstein condensation in certain 
cases. Later work of Toms showed how to interpret Bose-Einstein condensation 



in terms of symmetry breaking, in the manner of flat space-time calculations flll], |13 
The most recent study by Lee et al. [E9| showed how interacting scalar fields can 



be treated. Bose-Einstein condensation for self-interacting complex scalar fields was 
considered in ref. [pDfl . It is to be hoped that the techniques developed here can 
contribute to the development of this subject as well. 
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